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BOUNDARY-LAYER INSTABILITY AND TRANSITION 


SUMMARY 

The research conducted under this project has been directed toward the 
double objectives of providing (1) a rational foundation for the application 
of the linear stability theory of parallel shear flows to transition 
prediction and (2) an explicit method for performing the necessary 

calculations . 

The fundamental discovery upon which our subsequent work is based 
was that the solutions of the linearized, three-dimensional, incompres- 
sible Navier-Stokes equations u,p and the adjoint solutions u,p 
satisfy a "continuity" equation 


ll* 7 . J » 0 (1) 

where p is a pseudo-energy density (the dot product of u* and u) and 
^ is a pseudo-current. This result is derived and discussed in detail 
in Appendix A. 

We next considered (see Appendix B) the expansion of an arbitrary, 
two-dimensional solution of the linearized stream function equation in 
terms of the discrete and continuum eigenfunctions of the Orr-Sommerfeld 
equation in the half-space, y«(0,»D: that is, we considered boundary- layer, 

wake, jet or free-shear layer flows. Wc used equation (1) to derive a 
biorthogonality relation between the solutions of the linearized stream 
function equation and the solutions of the adjoint problem. Tiis is the 
biorthogonality relation for the mixed initial -boundary value problem. 


For the case of temporal stability, we used equation (1) to derive 
the formal solution of the initial value problem as a sum over the discrete 
inodes plus an integral over the contintium functions and showed that this 
expansion is conqplete. We found that the vorticity distribution at the 
initial time is sufficient information to determine the expansion coef> 
ficients and gave explicit formulas to calculate these coefficients. 

For the spatial stability problem, we showed that the continuum has 
four branches. We used equation (1) to derive the spatial biorthogonality 
relation and the formal solution to the boundary value problem. We have 
(see Appendix C) also derived the Fourier (in t), Laplace (in x) transform 
solution of the spatial stability problem and used it to show that our 
spatial expansion is complete. 

The boundary conditions for the spatial problem are the Fourier 
transforms, in time, of the stream function and its first three partial 
derivatives with respect to x, evaluated at x « 0. As it stands, this 
formal solution will not give a physically acceptable solution because, 
given an arbitrary variation with y and t at x ■ 0 of the stream function 
and its first three partial derivatives with respect to x, disturbances 
which lie on all four branches of the continuum will be excited. Therefore, 
as we show in i^endix B, the spatial wave packet will contain, in addition 
to waves propagating toward x ■ «, waves propagating upstream from x » • 
and standing waves whose amplitude increases towards x « <». 

A condition must be imposed that, for x > 0, all propagating dis* 
turbances are traveling in the positive x-direction and all standing waves 
have amplitudes which decay in the positive x>direction. It appears that 
this should be done by requiring that the stream function and its first 
three partial derivatives with respect to x, evaluated at x = 0, be 
orthogonal, using the spatial inner product, to all eigenfunctions on 
branches 2 and 4 of the continuous spectrum. 

It is easy to see that these two orthogonality conditions reduce the 
number of boundary conditions at x » 0 from four to two. This means that, 
for the spatial stability problem, the proper boundary conditions at x » 0 
are the specification of the temporal Fourier transforms of the velocity 
components u and v, for all y. Although these boundary conditions 
were derived from consideration of the continuum eigenfunctions, they 
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apply as well to the discrete, Tollmien>Schlicting inodes. We have not 
yet carried out a detailed investigation of the implications of imposing 
this orthogonality requirement on the boundary conditions: however, the 
immediate result that the boundary conditions at x • 0 are the specif i* 
cation of the temporal Fourier transforms of u and v for all y appears, 
on physical grounds, to be correct. 

We have presented preliminary numerical results of the application of 
this expansion method at the Fifteenth International Conference on 
Theoretical and Applied Mechanics (Appendix D). We considered the temporal 
stability problem and a simple initial disturbance. We assumed that at 
t • 0 the vorticity C was given by 

io.x 

C • Cq e «(y - y^,) (2) 

a periodic layer of vorticity at a distance y^ from the boundary. 

The stream function is then given by equation (55) of Appendix B, 
and it is easily seen that the expansion coefficients are [from equations 
(56a, b) of Appendix B] : 

' “o’ 

A^(a) . 6(a - 

The solution of this sinq>le problem, which is in effect the Greens function 
in y of the initial value problem, shows that the amplitudes of discrete, 
Tollmien-Schlicting inodes and the continuun functions are the products of 
the magnitudes of the corresponding adjoint functions, evaluated at y^, 
the height of the initial disturbance from the boundary and the vortex 
strengths. 

We applied this result to two different flows. The first is a slip 
flow past a bounding plane at y * 0. Although the base flow velocity 
does not vanish at the boundary, we required that the disturbance velocity 
vanish at y » 0. We found (Appendix B) that, because of the simple form 
of the base flow, all the calculations could be carried out analytically 


(3a) 

(3b) 


and the stream function could be expressed as a finite sum of exponentials 
and error functions. We found that the disturbance retains its identity 
as a periodic array of vortices for all time, but as time increases it 
diffuses, the vortex strength decays, and the centers of the vortices 
drift away from the boundary. 

The second flow we considered is the Blasius boundary layer. The 

velocity scale was taken to be the free-stream speed U and the length 

o 

vx/U . We chose a ■ 0.179 and R * 580.0. At this a and R, 
there are seven discrete Tollmien-Schlicting modes, one of which is 
unstable. We ‘numerically calculated the seven eigenfunctions and adjoint 
eigenfunctions and normalized them so that 


< <>_ 
n m 


nm 


(4) 


Plots of the amplitude and phase of the normalized eigenfunction and 
adjoint eigenfunction of the seven modes as a function of y, the dimen- 
sionless distance from the boundary, are given in Appendix D. These 
modes are numbered in order of increasing stability with mode 1 the 
unstable mode and mode 7 the most damped mode. 

The amplitude of a mode, say <|i^, excited by the vortex sheet at 
y » y^, is proportional to the amplitude of (>i^ evaluated at y^. It is 
clear from an examination of these figures that when the vortex layer at 
t a 0 is in the inner portion of the boundary layer, say y s 2.0 (the top 
of the boundary layer is at y = 5.02), there will be a relatively strong 
excitation of the discrete Tollmien-Schlicting waves. Modes 1, 2, and 
3 will have the largest amplitudes, and the higher modes will have sub- 
stantially smaller amplitudes. It is also quite clear that, when the 
initial disturbance i more than about four boundary- I aver thicknesses 
from the wall at t » 0, the discrete Tollmien-Schlicting modes excited 
by the disturbance will have extremely small amplitudes. We believe that 
this result is a theoretical explanation of the experimental observation 
of Kachanov, Kozlov, and Levchenko (1978) that vorticity disturbances 
passing above a boundary layer are very inefficient generators of Tollmien- 
Schlicting waves in the boundary layer. 


CONCLUSIONS 


We believe that we have created a rational foundation for the appli- 
cation of the linear stability theory of parallel shear flows to transition 
prediction and given an explicit method to carry out the necessary cal- 
culations. We have shown that these expansions are complete. We have 
also carried out some sample calculations which show that a typical 
boundary layer is very sensitive to vorticity disturbances in the inner 
boundary layer, near the critical layer; vorticity disturbances three or 
four boundary*layer thicknesses above the boundary are nearly uncoupled 
from the boundary layer in that the amplitudes of the discrete Tollmien- 
Schlicting waves are an extremely small fraction of the amplitude of 
the disturbance. 

After the completion of this grant we intend to continue these 
calculations. We will continue the calculations of temporal disturbances 
in typical boundary layers and begin calculation of spatial disturbances. 
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Expansions in Spatial or Temporal Eigeronodes of the 
Linearized Navier-Stokes Equation 

by Harold Salwen 

Department of Physics and Engineering Physics 
Stevens Institute of Tochnology, Hoboken, N. J. 07030 


The expansion of an arbitrary flow field in terms of the temporal or 
spatial eigenmodes of the linearized Navier*Stokes (LNS) equations for an 
incompressible fluid is developed from a unified perspective. It is i^hown that, 
for (v,p) a solution of the LNS equations for a given base flow and (u,q) a 
solution of the corresponding adjoint equations, a scalar “density", 

£{u»v), and a vec;^or "flux", ? (u,2f,v,p), may be defined such that^and r 
are bilinear in (u*,q*) and (v,p) and satisfy the "continuity" equation, 

3^3t + 7*r » 0. This equation is then used to derive biorthogonal ity relations 
between the eigenfunctions and adjoint eigenfunctions of the LNS equations for 
a general translationally-invariant problem. In the temporal case, the inner 
product is ■ ///»spu**v dr which is the natural extension of Schensted's 

inner product for two-dimensional disturbances and satisfies the requirements 
for an inner product in a Hilbert space. In the spatial case, the "inner 
product" is /// dydzdt which is not positive definite. The formal solution 
of the LNS equations is derived, in terms of the eigenfunctions and the initial 
or boundary conditions, for the temporal and spatial cases. It takes the form 
of the evolution of a three- or six-dimensional vector — (v , v , v ) in the 

A y £ 

temporal case or (v , v , v , 3v /ax, 3v /ax, p) in the spatial case. 

A y 4 y 4 


1. Introduction 


A few years ago, Grosch and I, after showing ‘that the Orr-Sonnerfeld 
equation for unbounded flows such as the Blaslus boundary layer possesses both 
temporal and spatial continuous spectra (Grosch and Salwen, 1975, 1978), set 
out to find the form of the wave-packet expansion for the temporal or spatial 
evolution of the stream function of an arbitrary two-dimensional "Irfinitesimal" 
disturbance In terms of the corresponding temporal or spatial el :'enf unctions. 

We sought to prove a blorthogonallty relation between the eigenfunctions 
of the Orr-Sommerfeld equation and of Its adjoint and, thereby, to solve for 
the coefficients of the expansion In terms of the Inner products of the adjoint 
eigenfunctions with the stream function at the Initial time or position. 

This worked out easily In the temporal case, with an Inner product equivalent 
to SchenstedS (1960) and only minor complications due to the Infinite domain 
and continuous spectrum. In the spatial case, on the other hand, we found that 
we didn't know the appropriate Inner product and v^e couldn't find any papers 
dealing with the problem. I therefore undertook the spatial expansion problem 
and, eventually, was rewarded with the result reported here— a unified treat- 
ment of the spatial and temporal expansion problems for solutions of the 
linearized Navler-Stokes (LNS) equations for an Imcompressible fluid. 

Section 2 Is devoted to the derivation of a “continuity” equation which 
Is used. In Section 3, In the definition of the Inner products and the derivation 
of the blorthogonallty relations. These, In turn, are used In Sections 4 and 5 
to derive the formal solutions of the (temporal) Initial value problem and 
the (spatial) boundary value problem, respecti vely. The application of these 
results to two-dimensional disturbances of a boundary layer has been presented 
In a separate paper (Salwen and Grosch, 1980). 

s 


In order for the formal solutions derived In Sections 4 and 5 to be 
actual solutions of Initial and boundary value problems, the eigenfunctions 
used In the expansions must form complete sets. Not all the eigenfunction 
sets one might want to use have ^ en proven to be complete but there are, 
by now, proofs of completeness for large classes of temporal eigenfunctions 
for bounded flows (Yudovich, 1965 and DIPrIma and Habetler, 1959) and temporal 
(Salwen and Grosch, 1980) and spatial (Salwen, Kelly, and Grosch, 1980) 
eigenfunctions for unbounded flows. 


2. •'Continuity'* eouatlon 

I start with the LNS equations for an Incompressible fluid with a base 
flow 0, 

v*v ■ 0 (la) 


3 V. -I 

3t’ * U-f V, ♦ V.f 1), J 


• Li 7^V. - ^ 


1 3X 


1 


, 1 • 1,2,3, 


and the corresponding adjoint equations , 


v*u ■ 0 


(lb) 


(2a) 


L ■ 5F ’ ’ * ^ 


30* i 


^ •Y 


3X. 


• u 


u , 1 - 1,2,3. 


(2b) 


* The complex conjugate, 0*, Is used here In order to obtain the correct formal 
expressions. In most applications, 0 will be real, so 0* « 0. 
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J. 




For any solutions (v,p) of (1) and (u.il) of (2), define* 


£ (u.v) ■ j ^ u* • V . 


3 ^ 

■ jr ^ ^ { *p |(tr u^*) • ^ ? •((! u^*)t 

» u [Ji,* r^v, - », 7^1, *j - |£- . j i 

with 

f (u.q.v.p) • (j pu**v) 0 + j 4 ^ [(V u^*) 

• u^*7 v^j ♦ j j^u*p + q* V , 


C(u»v) ♦ 7*f (u,q,v,p) ■ 0 . 


* The constant Factor, y s , Is Included In order to emphasize the 

e l 2 
and the energy density. ^ ^ v . 




* which has the form of a continuity equation relating the time Oerivatle of 
the "density", € , to the divergence of the "flux", r . For any fixed 
volume V bounded by a surface S , the "ccntinult; " equation (6) may be put 
Into Integral form, 

\ \ 

|r/// ^(u»v) dt ♦ // n«7 (u,2f,v,p) dS ■ 0 . (7) 

dt y 5 

3. ApollcatUn to a steady, translatlonallyinvarlant base flow . 

BlorthoQonalltv relations 

In this section. Equation (7) will be applied to the case In which the 
base flow and boundary conditions are Independent of x :nd t . For all x 
and t, the bate flow U(y, 2 ) , disturbance velocity and pressure lv,p) and 
adjoint velocity and pressure (u,q) are assumed to be defined In a closed, 
bounded area, A , of the y,z plane and to satisfy the boundary conditions 

v(x,y,2,t) ■ 0 , u(x,y,2,t) • 0 for (y,2)cC , (£a,b) 

on the boundary, C, of A. In this case, the temporal and spatial eigenfunctions 
discussed below will form discrete sets. (The extension to an unbounded area 
Is not too difficult (see, e.g., Salwen and Grosch, 1980) but It requires the 
relaxation of the boundary condition (8) and the consideration of continuum as 
well as discrete modes.) 

Because of the choice of base flow and boundary conditions, (1) end (2) 
are now Invariant with respect to translations In x and t and, therefore, 
possess solutions of the form 

7(x.y.z.t) • ifjt/.z).'*'** • (Ji) 


11 


(9b) 


p(x.y,z,t) a pQ(y,z)e^^®^ * 

u(x,y,z,t) « UQ(y,z)e^^^^ * (10a) 

q(x,y,z,t) » qQ(y.z)e^^®^ “ (10b) 

Because of (8), = 0 on S. Then evaluation of (7) over a thin slab 

perpendicular to the x-axis for functions of the form (9) and (10) gives 


i(v* - u») 


t(u,v) dydz • It II 


ax II ’’x 


• i(S* - o) II (u.q.v.p) dydz , (11) 

which will be used here to prove biorthogonality relations for the spatial 
and temporal eigenmodes. 


The temporal eigenfunctions are the solutions of (1) and (8a) having the 
form (9) with a real. These may be denoted by (^^p» P 2 ^n^’ corresponding to 
the x,t variation " “n(®)^). For each such solution, there is an adjoint 

eigenfunction (u^^ q^^) which is a solution of (2) and (8b) having the variation 
gi(ox - Vj^(a)t) v^(a) - o)*(o). Application of (11) to these functions 
gives 


(li)j^(a) - u^(a)) 


V 


am 




( 12 ) 


so that the integral vanishes when Ujjj (a) f appropriate 

normalization. 
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1 

4 


II ‘“'on.- iJ ■ 


03 ) 


(The expression is a function of x and t but is constant because the exponentials 
in the two factors cancel.) This biorthogonality relation for fixed a leads to 
the result for the full set of temporal eigenfunctions, 


<U 


am 


at 

, W = f II “sn> 




x=0 


dydz 




x=0 


dydz 


= 6(a-o) i5„„ for all t . 
mn 


(14) 


The spatial eigenfunctions and adjoint eigenfunctions are the solutions of 
(1), (2), and (8) having the forms (9) and (10) with u and v real. These may 


be denoted by C _ = (v „ P „)» with the variation and 

^ ain ujn , ujO 

^ m ~ (u „ <1 m)* with the variation As in the temporal case, 

vm ' vm, vm 


the eigenfunctions and adjoint eigenfunctions may be paired, with = aj^*(u) 

in this case. The analogous results to (12) and (13) are 


(a (lij) - a„(w)) [f (u „ 9 « V _ p „) dydz =* 0 

m n J J ^ ii)m , uin, ii)n, mn 

A 


(15) 


and 


1 r (u 0 V p ) dydz * 3 /2 tt 

•' X ' uiin.^mn, ain/mn' ^ inn'^ 


( 16 ) 
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which lead to the biorthogonality relation for the full set of spatial 
eigenfunctions 


03 



0 » 


• f f 

f 

^ q « V ^ p 
X w(ti» vn,'^vn 

J J 




dydz 


2tTiS(u)-v) 


(*u* q' V „ p „) 

X u)fn y uiHi y cjn y u)n 


t=0 


dydz 


5(u-v) for all X . 

mn 


(17) 


We thus see thaty with appropriate "inner products" <y-> and IC-I- 
the temporal and spatial eigenfunctions satisfy biorthogonality relations with 
the related adjoint functions. The temporal inner product <y> satisfies all 
the conditions ordinarily required of an inner product. The spatial inner 
product l-J y on the other handy is not positive definite. This is related 
to the fact that disturbances can propagate in both the downstream (+x) and 
upstream (-x) directions. 

4. Temporal expansion of an arbitrary solution of the L^S equations 
The temporal inner product introduced in (14), 


-u, v> 


u*(x,y,z) • V (x,y,z) dydz dx , 


(13) 
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is defined for any pair of ordinary vector functions of position. In 
particular, when applied to a solution (v, p) of (1), it involves the 
velocity, v, but not the pressure, p. It is natural then, in seeking an 
expansion solution in terms of the temporal eigenfunctions, to expand v 
alone in terms of the velocity part of the eigenfunctions. 

Let (v,p) be a solution of (1), satisfying the boundary conditions (8a). 
Assume that v can be expanded in the form 

00 00 

v(x,y, 2 ,t) * ^ I v^n(x,y,2,t) da . (19) 

n*l 

Then, by (14), the coefficients are 


c „(t) “ <u „ v> , 
an an , 


( 20 ) 


so that (using (7) ) 



1 (f It 



If 


v*r (u 


an, ^an 


Vf 


p) dydz dx * 0 


and 


( 21 ) 


c „(t) - c „(0) ® “^u \/> 

ixn an an , 


t = 0 . 


( 22 ) 
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The result is 


v(x,y,z,t) 



"“or. la 


(23) 


On the assumption that the expansion can be differentiated term-by-term. 


3X^ 


2 




+ V 






t=0 



3V _ 

r ””i 
L 3t 


+ U • vv 


an. 


+ V 


an 


7 U 




da 



(24) 


so that, except for an additive function oft only, 

oo 

p(x,y,z,t) = ^ 


OO 00 

f ^ 

<u,v>l p „(x,y,z,t) da . 
't=0 


n=l 


(25) 


Equations (23) and (25) are the formal solutions for v and p in 
terms of the initial velocity, v(x,y,z,0). 


5. Spatial expansion of an arbitrary solution 

The spatial inner product, [•I , introduced in (16), cannot be eval- 
uated in terms of the values of u , q , v , and p at a fixed x because 
involves x-derivatives associated \tith the second derivatives in (lb) and 
(2b). To get around this problem, one can regard the flow field at a given x 
as a 6-vector and make use of the fact that the velocities under consideration 
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A 


have vanishing divergence ((la) and (2a)). 


Let c and n be 6*vectors with components 


'2 • '3 ■ “z- '4 * • « 


5 *2 - 3x ’ '"6 


• p 


3U 3U 

= u^. = Uy. nj = u^. = 5/. ng « u^' * aT* '’e ' ^ 


and let 7 • u * 7 • v » 0 . Then, in terms of these components. 


Tj^(u,q,v,p) = % ^x 


|J*x*^3y * 3z”^ “ '^x ^3y * 3Z ^ 


3U * 3U * 

y . z \ 


+ (v U - u *v ' ) + (v U - U *V ') i 

y y y y ' ' z z z z 'J 


+ (Uj^*p + q*v^) 


so that 


00 

Jn, cjj * I II (u,q,v,p) dydz dt 


may be evaluated in terms of the components of i and n at fixed x . 
This choice of coordinates also eliminates second derivatives from (1), 
which becomes 


(30a) 


3 V 

— y. = V ' 
3x y 


(30b) 1 : 




(30c) 


av 

■ a V ' 
OX z 


3Vy‘ 

ax 


3U.. r 3v„ 3v 3v 3U 

+ |2.(^+U r-^+r-^v) 

u 3x X L n 'at y ay z az ay 7 ' 

a^v a^v. 




+ S. 


3U. 


V_ + ^ U V ' + — 

u az z y x y u ay 


(30d) 


!V 

ax 


0 >“z 


3U. 


y + £. ■ ■ y + 

u 3x X M 3y y 


r 3v, 
I (“ 

L li ^3t 


3V av au^ 

+ U r^+U ~+— 2 -v) 
y ay z az az 


aF“'J 




+ - U V ' + - 
u X z u az 


(30e) 


ax I ^'ay*- 


av 


a^v 

0^ 

< 

< 

ej 

< 

3U„ 1 

X\ 

- P 

(^+ U U ^ + 

'at y ay z az 

^v )1 

ax x'J 


Z5 

fT> 

1 

> 


, ■ y ) 

3y 

^ ^^x az " az '^z^ ” ^ 

(^ + 


az 


-) . (30f) 


It IS now straightforward to carry out the formal solution for the 
spatial expansion. The expansion is 


5(x,y,z,t) » 


2! I 

n*l “® 


(31) 


with coefficients 



(32) 
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Then 


dc * 

I f| 57''x (“un.ln/-'’) 


P 


57 '^y An 


kh ‘in.\n/'P> * It 8 'in/’l <« 


» 0 


(33) 


so 


c „(x) 

u)H 


C 




The solution is then 


(34) 


C(x,y, 2 ,t) = 


n»l 


IV 'I 


x=0 


^ _(x,y,z,t) dij 


(35) 
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ABSTRACT 


The expansion of an arbitrary two-dimensional solution of the 
linearized stream ^unction equation In terms of the discrete and continuum 
eigenfunctions of the Orr-Sommerfeld equation Is discussed for flows In 
the half-soace, y c C 0, -). A recent result of Salwen Is used to derive 
a blorthogonal Ity relation between the solution of tne linearized equation 
for the stream function and the solutions of the adjoint problem. 

For the case of temporal stability, the orthogonality relation 
obtained Is equivalent to that of Scnensted for bounded flows. This -elation- 
ship Is used to carry out tne formal solution of the Initial value problem 
for temporal stability. It Is found that the vorticity of tne disturbance at 
t * 0 Is the proper Initial condition for the temporal stability problem. 
Finally, It Is shown that the set consisting of the discrete elgenmodeb 
and continuum elgenfunctlcns Is complete. 

For the spatial stabUity problem. It Is snown that the continuous 
soectrum of the Orr-Sommerfeld equation contains four branches, "he b1 orthogon- 
ality relation Is used to derive the formal solution to tne boundary value 
problem of spatial stability. It Is snown that the boundary value problem 
of spatial stability requires tne stream 'unction and Its 'jrst three partial 
derivatives with respect to x be specified at x « 0 for all t. "o be apoll- 
caole to practical problems, this solution will require modification to 
eliminate disturbances originating at x • » and travelling upstream to 

X « 0. 

For the special case of a constant base flow, the method Is used to 
calculate the evolution In time of a particular Initial disturbance. 
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llipurbaacas to t bciadar/ liyar t.-.d, tharaby, iaitiatas trir.iltlca 


IS jsxacvr.. 


7ne General problem here is ehe solucion of ehe general initial 
and boundary value problems for disturbances to boundary layer flow— how, 
given .he fora of the disturbance at a tine t ■ 0, to find its variation 
with tiae and how, given the fera of the disturbance at all tines on a 
plane, x ■ 0, perpendicular to the boundar^*^ layer, to find out how it 
propagates dewnstrean. In this paper, we approach these problems, in the 
apprexiaation obtained by assuoing parallel flow and linearizing with 
respect to the disturbances, by eicpressing the solution as a sun over the 
discrete aomal nodes plus an integral over the continuum eigenfunctions 
of the Orr*Soraerfeld equation. If the (discrete plus continuum) 
eigenfunctions form a conplete set, this approach will yield a valid 
solution of the problem. 

Starting with H.aupt (1912), a number of authors have dealt with the 
completeness of the set of temporal eigenfunctions in a bounded domain. 
Haupc showed chat the eigenfunctions for two-dimensional disturbances to 
plane Couette flow torn a complete set and Schensted (1960) proved 
completeness for the eigenfunctions for two-dimensional disturbances :o 
plane Poiseuille flew and for axi-symmetric disturbances to Poiseuille 
flew in a circular pipe. Yudovich (1965) and DiPrima and Habecler (1P69) 
have proven completeness of the eigenmodes for a large class of bounded 
flows. We are unaware of any work on the completeness of the spatial 
eigenfunctions or, previous to this paper, on the completeness of the 
temporal eigenfunctions in an unbounded domain. 


2b 


In Section 2 we formulate the stability problem for two-dimensional 
disturbances to a parallel shear flow, U(y), 0 < *, in terms of the 

linearised equation for the stream function and boundary conditions. 

We next formulate the adjoint oroblem. A new result of Salwen (1979) is 
then used to derive a pseudocontinuity relation involving solutions of 
the linearised equation for the stream function and the adjoint solutions. 

This relation is then used to find tne general biorthogonality condition 
for wave-like disturbances to the flow. The bi orthogonal i ty relation is 
specialised to the cases of temporal and spatial stability. The orthogonality 
relation for the temporal stability problem is that deri*ed by Schensted (1960) 
and discussed by Reid (1965). 

The temporal stability problem is considered in detail in Section 3. 

The solution is Fourier anal^ced with respect to x. Then the formal solution 
of the initial value oroblem for the temporal stability of a two-dimensional 
disturbance to a parallel shear flow is expressed as an expansion in terms 
of the eigenfunctions. The expansion coefficients are determined y inner 
products between the initial disturbance and the eigenfunctions vf the adjoint 
equation. We show that the disturbance vorticity at t = 0 is the proner initial 
condition for the temporal stability problem. 

In Section J we examine the question of the completeness of the set of 
expansion functions for the temporal stability oroblem. Very recently, 
Sustavsson (1979) has treated the temporal initial value oroblem by using 
Fourier-Lapl ace transforms. He finds poles in the transform plane which 
correspond to the discrete T-S modes and a branch cut which corresponds to 


the continuous spectrum. We show in this section that the Fourier-Laplace 
transform solution of Gustavsson is ioentical to our Fourier transform, 
eigenfunction expansion solution for the initial value problem of temporal 
stability. We therefore conclude that our exoansion set is complete. 

The soatial stability problem is considered in detail in Section 5. 

The solution is Fourier analyzed in t. The formulae for the four branches 
of the continuous spectrum of tne spatial stability problem are derived and 
discussed. The formal solution of the boundary value problem for the spatial 
stability of a two-dimensional disturbance to a parallel shear flow is 
expressed as an expansion in terns of the spatial eigenfunctions. The 
expansion coefficients are determined by inner products between the boundary 
conditions at x = 0 and the eigenfunctions of the adjoint equation. The 
boundary conditions at x = 0 are discussed. We have not yet been able to 
prove completeness for the set of expansion functions of the spatial stability 
problem. 

In Section 6, we apply the results of Section 3 to the simple case of a 
constant base flow. In this case, we find the eigenfunctions and calculate 
and discuss the temporal evolution of a particular initial disturbance. 
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. The linearized, tvo dimensional Navier-Stokes eguaciotis 


Che 


biorthogonality relation . 

2.1 Foraulation of the problem 

The basic flow under consideration is a parallel shear flow, *J(y), in 
the semi- infinite region, y > 0. We are concerned with the temporal 
or spatial development of an '*infiritesimal”, rwo- dimensional distur- 
bance to this flow, (u(x, 7 , t) , v(x, y, t) , 0). In this case, u and 
V can be expressed in terms of a stream function, !?(x, y, c) , by 


u 



( 1 ) 


V 


3 # 

9x 


( 2 ) 


and the linearized Navier-Stokes equations reduce to a single partial 
differential equaticn. 
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In addicion 


ausc satisfy two boundary conditions at y ■ 0 




-v(x, 0, c) - 0 


( 5 ) 


and 



x.O.c 


u(x, 0, c) 


0 , 


( 6 ) 


nd a "fiaicaness" condition 
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As a consequence of eq. 


C). 




aust satisfy boundary conditions at 


infinity, 


3 ^ 

3x 



as y *• 


( 3 ) 


For fixed x and 


i?(x, y, t) belongs to a aanifold, M,ot functions , 


o(y), satisfying 
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:* 3 * 
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ail defined on lO, , 


(<5) 


continuous on lO, ») , 


( 10 ) 


3 . 


dy' 






3(0) • 0 and 3* (0) ■ 0 


( 11 ) 


and 


[ |3(y)|*iy and 

'0 


f" (dOi-. . 

rr— i 3y both exist. 

Jq 


(i:) 


The concinuun eisenfuncrions vhich will be discussed in Sections 3 end 5 
do not satisfy eq. (12). Instead they belong to a aanifold M* o M of 
functions satisfying eqs. (9) - ( 11 ) and a weakened condition, 

0 (y) and 7 ^ bounded in [ 0 , *). (13) 

cy 

Ve define an inner product. 


(f, g) = I f^(y) g(y)dy, (14) 

^0 

in M. The star denotes the coaplex conjugate. This inner product is 
defined for the full Hilbert space of functions satisfying eq. (12) and. 


.n that space, has the usual properties of inner products. 


The adjoinc problem 


?or functions g c M we define the adjoint, X , oz S in the 
usual way by 


(7^ 

.'jj 


{f (x,y. t) S"^g(x,y,t) ;dxdydt 


r; 


Ml w f(x,y,t)} {g(x,y,t) }dxdydt - Boundary Terms. (15) 

4 « J 


The definition of the adjoint ’used here yields an adjoint operator which 
is identical to the formal adjoint (Friedman, 1969, pp. 2,3). 

An adjoint stream function, 1? (x, y, t), is a solution of the 
adjoint equation (with U ■ U) , 


{Jj 

9^ — 



du 22# , 1 
dy 3x3y S ^ ^ 


0 . 
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wish Che boundar^^ conditions at y ■ 0, 



3vt 


x.O.t 


0 . 




and the finiteness ccndiffcn 





v|*}dy < 


(13) 


As above, equation (19) implies that ^ nust satisfy boundary conditions 




as 


y - 


(19) 


Vhen, as oelcw, ve look for solutions to the linearized stream 
function equation (3) which have a vavelike behavior in x and t* 
equation (3) reduces to the Orr-Soomerfeld equation and equation (16) 
reduces to the adjoint Orr-Somnerfeld equation. Our adjoint Orr- 
Sotmerfeld equation is the complex conjugate of the adjoint equation 
derived by Schensted (I960) and quoted by Reid (1965). The reason for 
this difference is that we define the inner product in the usual way, 
(14) while Schensted*3 definition of the inner product (f, g) involves 
f instead of f*. 
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.3 


31orchogonalicy 


Salv€n (1979) has shown chac cha solucions of the linearlocd, 
three ilaensional Navier~Scok.cs equations, 'u, p and the adjoint 
solucions u, p satisfy a "continuity" equation 
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3c 


J - 0. 
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where 
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-u)o4 : :(7uj)u^-u: c-v: 
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ind, IS before, :he scar denotes a coaplex conjugate. 

For the cvo dimensional disturbances considered here ve will 
introduce two new inner products. Let if be any solution of the original 
proolea and f be any solution of the adjoint problem, then define 
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'3x 3x 
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and 


S ! J dy . 


(24) 


with ch« X zoapontnc of J. Ualng (22) and axpresslag u, ? and 1, 
? in tens of $ and ^ , it can be shown chac 

- r 3“!^* 5^3iP 

• I "R"- 3 * Tx ^ ^ Tx T - 

^0 -ix 3x" 3x“ 3x^ 
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(2f) 


The fors of chase inner produces has been iecemined by the equations 
for the stream f'cnction and the adjoint stream function, however, we 
can use equations (23) and (25) to calculate inner products <f, g> and 
gZ , evaluated at fixed x and t, of any functions f(::, y, t) and 


s 

^ • f 


5(;t. y. :). 
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Ic is icrsighcfor<rsrd eo show chsc <f, g> is dsfinsd for :hs full 


Hilbsrt spscs of functions which satisfy aquation (*) and, in that spaca, 
has tha usual propartias of innar products. On tha othar hand, I f, gZ 
is not positiva iafinita. This is dua to t.ha fact that it is possibla to 
hava vavalika solutions to aquation (3) which propagata in aithar tha upst 
(•x) or downstraaa (-^x) diraction. 

Vith thasa dafinitions it is aasy co show that 

<!?,?' - -4- - 0 . C6) 

^ C « X 


for any solutions of tha original and adjoint problaas. 
If » and ara wava disturbancas of tha fon 
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aquation (26) raducas to 
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This aquation aay ba usad to dariva biothogcnality relations for tr.a 
aiganfunctions of both tha tatporil ind spatial stability proolans. 
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Tht ?«aporal 5c*blll:y ?ro;)l<a 


3.1 The Eigenvalues and Eiganfuncriona 

For :ht etoporal scabllicy problem we aodify ;he finiteness csndiiion, 
equation (T), to 


^'5x * ‘3y •• 


(’ ) 


This «nsurts zhAZ th« Fourl«r inc«gnl txpansion of ^ , 


’r' (.X. y, J ^^(y. :)e*^’‘ ia . 
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txijcs. If w« aisust that is of th« fonn 
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t.^en is i solution of the Oir-Sccner (eld equation 
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Slailarlv v« aaaumt chat cha adjoint solution, ^ . also taclafiaa 
aquation (7*) thua anauring that tha Touriar intagral txpanalin of 


(x, y, t) • j" V^(t, t)a^^* do 


axlata. It la aaauaad that ia of tha fora 


V",(y. t) - 5^(y)a 
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with tha solution of tha adjoint 3rr-Sceaarfald aquation 
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vith 


:* • -)*/o. 


Both and satiafy tha boundary/ conditions 


;^( 0 ) • ;^( 0 ) - 0 ^( 0 ) • 5 ^ 0 ) - 0 . 
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and aitha 
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if and arc In M, or che weaker condicion 

3* bounded as y ^ (43) 

^ U 3 Jl 

^ t 

if and arc in M . Those eigenfunctions which belong to M will be 

called discrete eigenfunctions. Those which belong to M* but not M will 

be called continuum eigenfunctions. 

It has been found (Mack, 1976, Grosch and Salwen, 1975, 1973a) chat, 

in general, there is a finite number of discrete eigenfunctions, \0 (y)} 

am 

with eigenvalues } and a set, {o ,}, of continuum eigenfunctions 
am aic 

with eigenvalues which depend continuously on a real parameter, k, 

in the range 10, ») . (Note that the k of this paper is equal to ak of 
Part 1.) 

The number of discrete nodes, which we shall denote by N(a), depends 

not only on a but also on R and on the form of U(y) and can, in some cases, 

be tero. The adjoint eigenfunctions also include a finite set, 

discrete eigenfunctions and a continuum, {3 . with eigenvalues ; and 

aK .«m 

respectively (see discussion following (30)). For a given k, 

^ • ikv 

and t , var^ like a linear combination of e * as y ». Ve there r ore 
ak 

find that 

'k-r: ^ 'k-^€ ^ 

0 , »(v)dk* and »(y>dk’ t M (4a) 

, ak ' j, , ak 
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and chac, for any square-incegrable f, 


^ do . u 

^ctk^ * ^ f^y ^ exist. 


(^5) 


Inner products between continuum functions, such as 

not exist in the ordinary sense but are definable in terms of the Dirac 
3-function (Lighthill, I960, pp. 10-21). 

The discrete eigenvalues must be searched for (Mack, 1S76), but the 
continuum eigenvalues follow from the asymptotic form - linear 

cotnbinacion or a““‘ ' ) of che eigenfunctions as y - » and U - ■ U(”) : 


(-k" - a")” - (iaRU, - iiLi.^j^) (-k" - 3 “) - 0 , (46a) 


'y 


(-k“ - 3")" + (ioRU^ - iSaj*j^)(-k” - 3“) - 0 , 


(iob) 


SO that both equations yield 


^3k * ^ 


(** 7 ) 


Ve also find that no continuum eigenvalue is also a discrete eigenvalue. 
Then 
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'"•’■A’ ^in’’ 


(43b) 


and 


3 0 , dk'> ■ 0 unless k, < k < k,. 

•^ak’ 5k' 


(43c) 


With proper Ubelliag and norraiizacion. i: is then possible to choose 
the eigenfunctions in such a way that 
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3.2 Expansion of an arbitrary ilscurbance. 


If the eigenfunctions fora a complete set, then, for any time, 

ve may expand ^ (y» t) as a linear combination, 
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N(Ct) 


(y, c) - I a (a. c)a (y) ! a (a, c)o, (y)dk 

a ^ , a an ] ^ •< ax 
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of those eigenfunctions . To find the coefficients <.a * and la, • we 

n xc 

make use of eq. (-*9) to cake inner products 
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We then find that 
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And, similarly, 
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where 

A (a) i a (a, 0) - *iK,(y. 0)> (5**a) 

a n an a 

A.^v.t) = ’ (5<ib) 

Then, referring zo equations (32) and (50), we find 

- r -t ' ► 
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and, similarly, 




(56 b) 


If Che discrete and continuum eisenfunctions form a complete set, 
then equation (53^ constitutes an expansion of the stream function of an 
arbitrary disturbance in terms of the discrete ^Tollmien-Schiichting) 
and continuum wave solutions, 






ok 


of :ae iiscurbanca equacioa, (3). vi:h coefficients iecernined by the 
initial fora of the disturbance ^x, y, 0). In the next section «e will 
Show that the discrete and continuum eigenfunctions are a complete set. 
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4. Completeness of the Temporal Expansion Functions 

Gustavsscn (1979) has carried out a formal solution of the initial value 
problem of temporal stability for three dimensional disturbances. He uses 
the same coordinate system as we do with the addition of the z coordinate 
in the cross stream direction. The formal solutionis obtained by taking 
Fourier transforms in both x and z and a Laplace transform in t , formally 
solving the Orr-Sommerfeld equation in the transform space, and formally 
inverting the transforms. If we eliminate the z-variation of Gustavsson's 
solution and his Fourier transform in z (replacing his k by |a|) the two 
solutions should be identical. Both Gustavsson and we express the solution 
in physical space as an inverse Four.ier transform over s, the transform 
variable in the x direction. In order to show that these two methods yield 
identical results it is therefore necessary to show that his formal solution 
in Fourier space, v as given in (G13) , is equal to the factor in curly 
brackets in our equation (55). 

In order to do this we must first translate Gustavsson's notation into 
our notation. Setting 2*0, after (G3) it is easily seen that we have the 
following correspondence, 



baper. 


This paper 

Gustavsson 


i v/a 

|a| 

k 

Ui 

is 


k a 

U, 1 

1 n (G3) and thereafter. 

Gustavsson gives the formal solution in Fourier space in equation (G18). 

It consists of a sum of the residue values at the poles plus a contour 
integral along a branch cut. Using the definitionsof W as the Wronskian^ 

Aa4 the D-, given after (Go). and the 5 ., equation (G7), it is quite straight- 
forward to show that the residue, at a pole s^ is 

= (e^vVw) lim i'(s-s^) Ca^ (s)o^ (y.s) + a 2 (s)o 2 (y .s)]). (58) 

^ ^ V 

Therefore the residue consists of a linear combination of >•] and ^2 
the solutions of the Orr-Sommerfeld equation that approach zero as y - 
i.e. they satisfy (34) and (42). At s » s^^ , *1 and >2 satisfy the usual 
eigenvalue condition at y * 0 for the discrete modes of the Orr-Sommerfeld 
equation, condition (Gb) (at the bottom of oage 1503). This linear comoination 
thus satisfies (41). Therefore the residue at s is prooortional to our discrete 
eioenfunction : (y) with eigenvalue m , and 

XJ ^ IV 


48 


It is well known (Coddlngton and Levinson, 1955, p. 101, problem 19) 

★ 

that COj/W] , the complex conjugates of the functions used in (Go), are 
solutions of the adjoint equation (39). It can be seen from the form of 
(Gil) and the definition of our inner product (23) that 



j • 3,i, 


(59) 


so that tne coefficient of 5 ^^(y) in the residue is the i.nner oroduct of 
some solution of the adjoint equation with ;(y,0). Finally, some straight- 
forward, but tedious, algebra shows that the particular linear combination 

♦ 

of the 0. involved satisfies the boundary conditions (41) and (42) and 

J 

therefore is a multiple of our s . We thus find that tne residue at s is 

av V 

R * d A (a) ) e , (60) 

V av V av 

with d^^ independent of y and A^(a) given by (54a). Before determining d.^,., 
we turn to the contribution of the branch cut. 


Using the ^act that our o» * is, it is clear from (G14) that the branch 

cut in the complex s plane is our continuous soectrum in the complex 

a) (or c) plane and that the brancn point, * 0, corresconds to the limit 

2 

point of our continuous soectrum at c » U.| - i a /R, with » 1. The function 
f(a,k; y) in (318) is, by (317)and (G19) a linear combination of the solutions 
of the Orr-Sommerfeld eouation which are, as y - asymototic to 
e'’’'-'^, and It can be snown, using (G19), (320), (321 ), and (322), that 

F(a.k;0) » (^) * 0 , 

°^y«0 


( 61 ) 


and so F(i,k; y) is some multiole of our continuum eigenfunction Sj|^(y)* 
Further, It Is obvious that, In (55), 




( 62 ) 


In (G18) with » 1 . 

Just as for the discrete modes, the {ah , v » 2,3,4 in (G21) are tne 

'J 

inner product of some solutions of the adjoint equation witn Using the 

definition of the E._ in (G22) and the definitions of the ID_ as given 

in tne next to last paragraph on oage 1604^ some algebra shows that the 

particular linear combination satisfies the boundary conditions at y » 0 

and so the inner product in (G20) is a multiple of the inner product of our 

continuum adjoint, hk^^^’ initial condition. Therefore, the integral 

term in (G18) is 
*• 

^ ' *^ik dk , (63) 

b 

with A|^(a) given by (54b) and d^j^ indeoendent of y. Gustavsson's result (G18) 
thus takes the form (in our notation) 


(y.t) 


N(x) 


d A (,) ; (y) e'^'^-.v^ 


b 

Both Gustavsson and we may cnoose our initial condition arbitrarily, 
provided that the various integrals of this function with the adjoint functions 
exist. I* we cnoose the initial condition tnat i»^(y.O) is one of 
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the discrete e1 jenfunctlons , say *'j^(y)t then In (55) 


• -m 


A|^(a) • 0 


(65a) 

(65b) 


In Gustavsscn's formulation (63) we have 


A (a) - 5 „ 

3V V VlH 


A|^(s) 5 0 . 


(66a) 

(66b) 


We thus see that 


d ; 1 

IV 


(67) 


If we then choose the Initial condition 


k+e 

f 

' k-e 


5 i,i(y) dk' . 


( 68 ) 


a similar argument shows that 


^k ^ ’ 


(69) 


Substitution of d • d . * 1, (67) and (69), makes eq. (64), derived from 
Gustavsson's solution. Identical to the curly bracket in our exoanslon 
solution (55). We have tnus snown tnat the formal solution obtained oy 
Gustavsson frc.n t ,e Fourler-iaolace transform is identical, term by term, to 
our formal expansion solution. 
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since doy souare-integrable solution possesses a Fourler-Laplace 
expansion, we have shown that our expansion (55) Is complete whenever It 
Is valid to separate the Fourler-Laplace transform solution Into a sum over 
the poles plus an Integral, over the branch cut— that Is, whenever the sum 
over the poles (discrete eigenvalues) converges. This Is, of course, Oso 
the condition for the validity of Sustavsson's solution. 

For the Blaslus boundary layer, the numerical evidence (^lack 197f) 
Indicates that, at a given R and a , the number of discrete modes is finite , 
so that the sum over the poles Is a finite sum. If this is so, then tne 
above condition is certainly satisfied and our expansion functions form a 
complete set. 

We have shown that the Fourler-Laplace transform result and the eigen- 
function expansion result are different forms of the same solution of the 
Initial value problem to be chosen according to convenience In a particular 
case. The eigenfunction expansion formulation gives explicit formulae 
(54 a,b) to calculate tne expansion coefficients. This allows one to calculate 
the amplitudes of tne discrete modes (TS modes) and the continuum functions, 
given the Initial distribution of vorticity. 
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5 


Tht icabtllev ?rcbl«a 


5.1 Th« «ig«nvalj«s and aiganfunceions 

Tha fialc*nasa condisicn, aquation (?), la aodlfiad ior tha tpaclal 
acabilltv problaa to 


JO 


Jr 


3y* 


- -^r.-dcdy < 


(:••) 


rhla asauraa that tha Fourlar Intagral axpacaion of ^ 


(x, y, c) ■ (x, y)a"^‘*'^du 

J -a ** 


(70) 


exists. If vt issuce thsc ^ Is of tne fora 


H^Jx. y) - 


. . » i 2 x 
: (y)a 


(71) 


then : is the solution of the Irr-Scnnerf eld tqustlon 


{1‘ - lilC(iL’ - J)1 - -J. ^:*t - D 

5 - iy* - 


(72) 


with givan by (16). 


Si:nilarly, we assume chat Che adjoint solucion, ^ , also satisfies 
equation (7**) thus ensuring chat 




( 73 ) 


exists. Ve assume that 


y) 


0^(7) e 




(74) 


Then 0 is che solucion of Che adjoint Orr-SoTnmerfeld equation 


LL^.* «[(=*» - »)L,. + f i:>;^ 


(75) 


The boundary conditions are 


5 (0) - o'(0) - 9 (0) - j'(0) - 0 , 

Jj ' U) JJ (JLt 


(76) 


and 


j-j'-J-a’-O as y-*. (77) 

Jj ^ ^ 


if p and P are in M, or 

jj 
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: . 


bounded is 


y ® » 


( 78 ) 




^ I 

jj 


. ^ » 

ir b and b are in M . As above, che eigenfunctions which belong to 

jl) uJ 

M are the discrete eigenfunctions and those that belong to M but not M 
are the continuum eigenfunctions. 

jordinson (19*^1), Corner, Houston, and R.oss (1976^, and Murdock and 

Stewartson (1977) have shown that there is only a finite set of discrete 

eigenfunctions, {o (v)}, with eigenvalues {ot }. The set of discrete 
^n jjn 

adjoint eigenfunctions, [o }, with eigenvalues (a ; is also finite. 

ojn ^n 

The number of discrete nodes, N(4j) depends on R as well as j and can be cero. 

In part 1 we shoved that, in an smbounded domain, the spatial 
stability problem always has a continuous spectrum. Since then we have 
discovered (Grcsch and Salwen, 197Sb) , that the spatial continuum of 
Part 1 is only one branch of a four branched spatial continuum. It is 
quite easy to show the existence of the four branches of the spatial 
continuum. Ve look for solutions to equations (^2) and 5 and 

ojiC 

^ ^ikv 

O^^(y), for a given real k, which vary like e" * as y » (the k used 
in discussing the spatial continuum in Part 1 is 2. R times the k used here). 
Noting that, as y ■*• », U U^, a constant, and U’ , U 0, we have 

(-a" - - k* - ^ i^R) - 0 . (79a) 

and 

k”M-i*^- k' - ia*RU^-i-R) • 0. (79b) 

(Note that equations (^9a and 2 ) are complex ccn jugatas . 
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J 


is obvious cha: there are four roots, 1’*^;, j 


and 1, the roots of 


oi“ * - k’ - i«R - 0 


- ilc , a, - -ik 


The eigenvalue the root zi equation (,80; vich positive real part, 
is the continuua eigenvalue discussed in Part 1. As was discussed in 
Part 1, the eigenfunctions of this branch of the spatial continuun 


are waves 


propagating in the downstrean (-^x) direction and decaying in 


amplitude as tney travel. In the same way it can be shewn that a,, is 
the eigenvalue of a continuum eigenfunction which is a wave traveling in 


the upstream (-x) direction and decayir 


ying as it travels, 


The tree stream speed, can be taken to be unity for a bcundap.' 
.ayer, wake, or tree shear flow. In most cases of interest j/R < < 1. 
It is easy to see that, with U, - 1, and j/R < < 1, 


1 - - '3 

a, , • t(oj/v> . 1RC4(1 • Y) - .vj/R)* 




J 


■< l/'> 

V - (1 + 4k' ^ . 


Define, as usual, the phase speed by 


★ * ; I 

c, - a j/ a,' 

J J J 


Then as k 0, 


. ^ 2 , 2 - 
j iRL (iA)* ^ k")/R*-* 




a, - -a) - iRCl ~ k”)/R'I, 


' 1 - t[(uj' * k')/R'I/U/R) , 


: -(uj/R)' i(WR)'l + * k'), R'; 


%"hile as k - *, 


; (^^R/2k) ik , 


1, > . 


oj'R/Tic'’ - i-»/k . 


(83) 


(34) 


(35a) 


(85b) 


(86a) 


(36b) 


•.37a) 


(87b) 


(38a) 


(88°) 


The damping rate, for Che spatial eigenfunctions, is Im(ct) and the 
phase speed is ReU). rhe equations given above show that the eigenfunctions 
on branch of the spatial continuum, for boundar^/ layers, wakes, and 
free shear flows, always have both a very large damping rate and a very 
small phase speed. This is in marked contrast to those of branch one, 
which, as was shown in Part 1, or can be seen from the above results, 
contains lightly damped eigenfunctions some of which have a very slow 
phase speed and seme of which have a phase speed nearly equal to the free 
stream speed. 

The spatial continuum eigenfunctions of branches 3 and 4 are standing 
waves in x because they var;.* like 


ictqx -k:c 
e - e 


l:iAX 


e 


(89a,b) 


As In she cemporal case, the Inner prcducss between the spatial 
continuum eigenfunctions do not exist in the ordinary sense but can be 
defined as 5 functions. Then^with proper labeling and nomalitation, 
it is possible to choose the eigenfunctions such that (vita the super- 
script i or J indicating the branch of the continuum) 



3 • _ 
un 



(90a) 





(90b) 





(90c) 
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; * ^ * -i 

jjn -un 


-t(a » a ) 

'{C 

„ R -in 

0 


do do » 

’‘^n “ '^-jn ^-n' “ iy dy 


+ uC(a” • + i . s * a i ' 

jti -in -in -;n -jn 


do do » 

^ ^ : -:}iy . 

dy dy dy -in 


and there are analogous expressions for the 


inner oroducts in (90b) 


and (90c). 


j 


5.2 Expansion of an Arbitrary Disturbance. 

If Che spatial eigenfunctions fora a complete set, then, for any 


x)5^*'(y)dk. (92) 


(X. 

Jj 

y) as 


<4 


a,(-), x) ; (v) - ; 

a*l 

^ i.i 


In order co find the coefficients'ia (^, x)} and (j, x); ve use 

n tc 

equacion (90) to sake the inner produces 


N(uj) 

: t (x. y)3 ■ - a I (a), x)o ' - a (u). x) . 

-jn j * , n nn n 

n -1 


(93a) 


V' (x. y)2 - I a,^-l^(^. x)o(k - k')d... dk' - x) . (93b) 

^ J Q iC •J 


Then 


3a (a), x) 
n 

3x 


^n 


0^ 
jj- 
— . 
5X 


ia :7 , : 

u-n -jn uj 


in a (u), x) . 
ujn n 


(9da) 


and 


oO 


# 


(u), x) 
3x 


3'^ 

* ^ 

“ 3x " 


^ ^k’ ^y 




(94b) 


30 chat 


a (uj, x) 
n 


V 

A (uj)e -*'' . 

n 


(95a) 



(u). 


x) 


,(i) 


(oi) a 


’4k^x 


(95b) 


where 


A^Oo.) S a^U, 3) - I • 


(96a) 


*.(«(.) = 4‘>ta, 0) - ;5<|;', V'/o, y)-. 

fC ** 


(96b) 


Fron equations (73). (92). and (95). ve have the fental solution to 
the soatial stability oroblea for the r.o dinensicnal. linearited 
Navier-Stokes equations 


JUa) 


f (X. y. t) - " { : 


ia^n-'^ 


. -t.t . 


i-1 ^3 


“x "■“ '-k 


(97) 


bl 


1 


Define 




1 ^( 0 ) 


(y) 


5 ^'iO, y) • ! $'(0. y. t)e^‘^'dt . 


(98a) 


,,v (x, y) , . , _ . 

^ <-^> • * ! t ’^^V , ‘--- (58b) 

x»0 ' -» x»0 


,3$ix. V. c)^ i 


uJC 


>^<=>(;0 = (. 

SjJ 


^> . . J, . 3-fe. V. .1 , 

x *0 :c"0 


iijJt .. 

e dc. ( 98 c) 


U/( 2 ) 


(y) i 


5 '^ I/' Cx, y) 

( ) 


^ ( 


sx 


x *0 


J^!?(x. V, c) , 
, 3 

3X 


e"''' dc. 


x -0 


( 98 d) 


Then 


A (j) - : 5 'f' ( 0 . y)Z 

n jj 


••• V T 

Jo ^ -n J 




. J* 

wjn ^ 


C) 


, : ;. y.i) 

-n «n ^ j 


iz ^ r - i:n -7^ ^ i 

“ “ dy :y 


-jn ^n jj 


d 3 * df^ 

JIl J 


1 ) 


iv IV 


' ’ > 

**n 


Li>( 0 ) . 




'^v 2 ) _ > 4 ;U' _ 2 ‘ 0 ) 


di ^ 
iv 




^ ' 3 ':;. iy. 


av 


( 99 ) 


t.: 


i 




( 1 ) 

and ch«r€ is a iiailar axprtssion for (dj) • 

This is Che foraal solution of the spatial stability probltni 
for an arbitrarily iaposed boundar/ condition at x ■ 0. The boundary 
conditions which aust be specified are the Fourier cransforas in tiae, 
of the stream function and its first three partial derivatives wish 
respect to x, evaluated at x - 0. 

As it stands, this formal solution will not give a physically 
acceptable solution because, given an arbitrary 5Po, y, t) and derivatives, 
disturbances which lie on all four branches of the continuum will be 
excited. Therefore the solution will contain, in addition to the waves 
propagating towards x ■ * and the standing waves whose amplitude decays 
towards x - waves propagating upstream from x • « and standing waves 
whose amplitude increases towards x ■ *. 

A condition must be imposed that, for x ^ 0, all propagating dis- 
turbinces are traveling in the positive x direction and all standing 
waves have amplitudes which decay in the positive x direction. It appears that this 
should be done by requiring that ‘r' (0, y, t) and its first three partial 
derivatives -with respect to x be orthogonal to all eigenfunctions on 
branches Z and -4 of the continuous spectrum but we have not yet investigated 
the implications of imposing this condition on the disturbance stream 
functicn at x - 0. 
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6. Applliitlon CO :h« 7«niporal Dtvelcpmt nc at a Modtl "low 

In this section, ws apply ch# results of section 3. to the siaple 
base flow, 


L'(y) ■ ■ constant, 


' 0 , 


( 100 ) 


which is a slip flow pest a bounding plane at y ■ 0. Though the base 
flow velocity does not vanish at the boundary, we still require the 
disturbance velocity to be zero at y - 0. 3ecause of the siaplicity or 
the base flew, the expansion functions are elenencary functions. 

In 6.1, we find the expansion functions. In this case, there are no 
discrete eigenmodes; all of the eigenfunctions are continuum functions. 

In Subsection 6.2. w« »ol'^« ti=e-d«veiopcent of a particular 

initial disturbance oy expanding in terns of these eiger.f’unctior.s . .he 
initial disturbance chosen is a periodic layer o. .or.-»--/ w-r..»ned 
plane parallel to the (y - 0) boundar.'. Because of the staple tom of the 
initial disturbance and the siaplicity of the base .*ow, .. is pcssib.e 
obtain the solution in closed fora in teres of error functions. 

6.1 The E jer. functions 

-or the base flow of Equation(lOO) the .fferential equation. f30. 
for the expansion functions becoaes 


(-i-r - n' - in.^ (V, - t)U-=-T - -■ 


3 . 


( 101 ) 


dy 


i 


vlth :h« gtneral solution (for i f 0, c • U^), 


At 


-la y 


^ 3t' 


Ct 


?y 


D« 


-?y 


( 102 ) 


uhtrt 


2 

? 


a"* 


laR(U, 


c). 


( 103 ) 


(In :hi* cas« of conscanc 'J, 5 aujc satisfy th« saat iifJerenclal tquacicn) . 
•n addition, i aust satisfy iqs. (11) and (13); i,a«, 5 and 5 aust 
vanish at y - 0 and b« boundtd in !3, ») . Since t ^ ^ is unbounded, 

3-0. To satisfy t.he boundary condition at the origin, we aust then have 


“ 1|V . 2 

Ate • - cosh ?y ♦ sinh ?y. 


( 104 ) 


which is unbounded as y • unless p is purely iaaginar^?. The solutions 
are then given by 


? - ik; 0 < k < » , 


( 105 ) 


"ik 





3K 


(y) 




cos !cy) 



< 7 .. 


( 106 ) 


I h5 

1 


» 


vh«rt :ht 




.izaiion conacanc, 


A 


:ik 




no7) 


is dscaralnsd by :h« csndicisn 

<a , , b . - 5(k - (108) 

3lc bk 

In this case, where the and are known explicitly, one -nay 
show directly that, for F(y) any continuous, d1 f'^erentiable , square- 
integrable function In [0,*), 

m 

' *^ik' (109) 

0 

thus confirming that the set of complete for functions In 

with P(0) • 0. 


bo 


6.2 The Tecpcril Ivcluclsn of An Inicial Diacurbanct 

In Criaar :o itaonacraca :ha applicacion of ciUa axpanaion ctchnlqua, 
va eonaiiar c.ba particular iaiclal iiacurbanca 


lP(x, 7, 0) - T(x, V, 0) - Tq i(y - 7g). 


a partodic layar of vorticity at a diacanca y^ froa tha boundary. 
Following aaction 3.1, wa find t.^ac t.ba scraaa f-jnetion ac any tiaa will 
ba givan by 


(x, y, t) 


. ^ J 0 * 


where 


A.^(3) • , T(». F. ixdy. 

*■ * ' • 0 ' -* 


Ii ij ea5il7 seen, by sucstitucing Eq. (HO) ints Eq. (11?) 




to :hei 




^(x. y, c) - *'^'‘ 


da 


^iao(x - Uic) ^-3^c/R 


Q 


Aicer using Sqs. (106) and (107), for ^ we find chat each tern in che 
integral is expressible as sums of error functions. The results are given 
in an Appendix. From these results « it can be shown that, for c - » with 
7 fixed, 




e 

/ t 


cos :^x - Uj^t) • (function of y) 


(115) 


and, for y * vith c fixed. 


cos :^(x - Uj^c) • (function of t) . 


(116) 


It is clear that, even though the indi*/idual eigenfunctions used in the 
expansion oscillate with constant aaplitude as y * *, the wave packet behaves 
like e as y -* *. 
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Figure 1 shows contour plots of the stream function for the disturbance, 
in a frame of reference moving with the free stream velocity, at six different 
times. Ve have chosen ■ 1.0 and y^ *1.0 for the example shown here. 
Contours of the disturbance stream function have also been calculated 
for other combinations of values of and y^ and, for these other values, 
the evolution of the disturbance in time is quite similar to that shown 
in figure 1. 

In figure 1 the (-^) and (-) indicate the position of the maximum and 
minimum values of the scream function. These maximum and minimum values 
are given in the caption to the figure. The flow is counter-clocicvise 
around a maximum and clockwise around a minimum (-) . 

It is clear from this figure that the disturbance, which is a periodic 
vortex sheet at t ■ 0, retains its identity as a periodic array for ail 
time, but as time increases it diffuses, the strength decays, and the 
centers of the vortices drift away from the boundary at y ■ 0. 

Ve could, of course, generalize this model problem by considering an 
initial vorticcy distribution in the y direction. ^*e have not carried out 
this calculation because cur intent in solving this model problem was to 
illustrate the expansion procedure and ve do net think that it %’arrants 
further elaboration. 
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Figure L 
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Contours of the disturbance scraaa function for the nodel 
problem in a frame of reference moving with the free stream 
velocity at six different tir4es. In this example * 1.0, 

^0 * ^0 * “^ere are tventy contour lines on 

each plot. The values of ^ on these contours are 0.95 i^max, 
A max, ... , *0.95 max. The and ( — ) indicate 
the positions vhere ^ 

. (1/ 




O . 

• min 


max 


(a) :/R - 10*^ , ^ . 0.-25 

“^max 


(b) t/R - 10~“ , - 0.359 

“ aax 

(c) t/'R - 10~^ , ^ - 0.212 


aax 


(d) c.R - 1.0 , » 0.223 :< lO”^ 

max 


(e) t. R - 5.0 , ip, 


niax 


O.IOS X 10 


(f) C/R - 10.0 , ^ » 0.3S3 :< iO"’ 

r.ax 


'^1 


•0 



Appendix; Solution of the Model Problem 

In section 6 we showed chat Che scream ^unction for the model 
problem is, equation (114), 



y. c) 


• ^1^) ^-ciqc/R 


7 

•k^c/R 




vhere (y) is 

0 

T - C/X, 


given by (106) 2 nd (107). Substituting for 2 nd 

s 

integral it is straightforward 0 shew that, with 


y. t) - t) - e^O^ r,(a^, t. y^) 


+ a„ 13(03, T. yg) - e "^^0 1 ,( 5 ^. T, y) 


2' O’ 


+ i 1,(3-. T, y + y-) + i I,(ct . T, y - y-) 


2g 13(3^. T. y - yg) - 3g e-'O^O 13(^3, r. /) 


0' 0 


2 '‘0 ' 


- ^ 0 ^ -2 -0 •’ ^ 


(A 


• f 
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here th( 


i functions Ij are given by 


2 -a 


■ ? • 


■t r 

J n 


^ 2 

Oc’ + i) 


-3dk 




, ^ 2 
■^) - ^ ’1 . 




(A2) 


_ , - 2 -a^T 

I, (a, 2) = - e 


r" .2 ,2 

CCS k-dk 

*'o (’< + a")- 


C(1 + 2a‘T - a2)e"'^^ arfc (ax'* - 


(1 


-h 2oi^T + a2)e^^ erfc (ai^ aZ/ZctT*) 


U 

-a"T -a"2“/4a"T^ 
€ e j 


(A3) 




sin \<Zdk 


- CCa^T + :tZ)e^'' 


erfc 




- Ca‘“T - a2)e*’"“ erfc (::T ‘ - aZ/ZciT-)-. 


(A4) 


1 


1 


cos kZdk 




-< T 


2 2 Z'* 

(k^ + r)" 


• ^ C(1 - :a^T + aZ)e"'^^ erfc (ar’* - 32/2 citS 


^ {I - :a*"r - :sZ)e’"^ erfc (ax * aZ/ZaT*) . 


^ -» -) 2 "* 

^ oax * -a"x -a*Z /4a"x-, 
+ e e J 


(.^) 


and, as usual 


erfc(Z) 


r-.-r 


r'T * : 


ac;. 
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ABSTRACT 

In this paper we give the solution of the boundary- layer receptivity problem: 
that of determining the amplitudes of the Tolimien-Schlichling modes and 
continuum eigenfunctions of a boundary layer given the form of the velocity 
profile and the disturbance, within the context of incompressible, linear 
stability theory for a parallel shear flow. We give the formal solution 
to the initial value problem for temporal stability and give the proper 
initial condition for this problem. The formal solution of the spatial 
stability problem is also given and the proper boundary conditions at x • 0 
and radiation conoitions at x ■ are discussed. We give examples of the 
application of this method to the calculation of the temporal evolution of 
a particular disturbance in two flows, a constant base flow and the Blasius 
boundary layer. 
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and 
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Department of Physics and Engineering Physics 
Stevens Institute of Technology 
Hoboken, N.J. 07030 


SUMMARY 


The last ten years has seen an increasing use of the theory of hydro- 
dynamic stability to predict transition in boundary layers. Mack (1977) 
gives an excellent, up to date review of various transition prediction 
methods. All of these methods include at least one unknown parameter 
Aq, the initial amplitude of the disturbance in the boundary layer. 

There are numerous discussions of the boundary- layer receptivity oroblem, 
that is, the problem of determining Aq given the velocity profile of the 
boundary layer and the disturbance (Obremski, Morkovin, and Landahl, 1969; 
Mack, 1977; Berger and Aroesty, 1977). All of these authors conclude 
that the mechanism by which free-stream vorticity and sound disturbances 
generate Tollmien-Schlichting waves in a boundary layer is unknown. 

In this paper we give the solution of the boundary- layer recep- 
tivity problem within the context of incompressible, linear stability 
theory for a parallel shear flow. The expansion of an arbitrary two- 
dimensional solution of the linearized stream function equation in terms 
of the discrete and continuum eigenfunctions of the Orr-Sommerfeld 
equation is discussed for flows in the half-space, y e[0, ») . A recent 
result of Salwen is used to derive a biorthogonality relation between the 
solution of the linearized equation for the stream function and the 
solution of the adjoint problem. 

For the case of temporal stability, the orthogonally relation 
obtained is equivalent to that of Schensted (1960) for bounded flows. 

This relationship is "sod to carry out the formal solution of the 
initial value problem for temporal stability. It is shown that the 
vorticity of the disturbance at t = 0 is the proper initial condition 
for the temporal stability problem. 

For the spatial stability problem it is shown that the continuous 
spectrum of the Orr-Sommcrfeld equation contains four branches. The 
modes on these brances are (1) waves propagating doi>mstream, (2) waves 
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propagating upstream, (3) standing waves whose amplitudes decrease down- 
stream, and (4) standing waves whose amplitudes decrease in the upstream 
direction. The biorthogonality relation is used to derive the formal 
solution to the boundary value problem of spatial stability. We show 
that the boundary value problem of spatial stability requires the stream 
function and its first three partial derivatives with respect to x be 
specified at x = 0 for all time. The imposition of a radiation condition 
downstream, i.e. at x = ®, eliminates disturbances which originate at 
X = « and travel upstream to x = 0. The imposition of this radiation 
condition reduces the number of independent boundary conditions at x = 0 
from four to two. 

We give two examples of the application of this method to calculate 
the temporal receptivity of boundary layers to a disturbance. We specify 
the disturbance at t = 0 to be a vortex sheet parallel to the boundary and 
sinusoidal in the streamwise direction. We then calculate the evolution 
in time of this disturbance in (1) a constant base flow, for which the 
calculation can be carried out analytically and (2) in the Blasius boundary 
layer for which we calculate the amplitudes of the discrete Tollmien- 
Schlichting waves and of the continuum eigenfunctions numerically. 
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